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To formulate the zero modes in a finite-size system with spontaneous breakdown of symmetries in 
quantum field theory is not trivial, for in the naive Bogoliubov theory, one encounters difficulties such 
as phase diffusion, the absence of a definite criterion for determining the ground state, and infrared 
divergences. A new interacting zero mode formulation that has been proposed for systems with a 
single zero mode to avoid these difficulties is extended to general systems with multiple zero modes. 

It naturally and definitely gives the interactions among the quantized zero modes, the consequences 
of which can be observed experimentally. In this paper, as a typical example, we consider an atomic 
Bose-Einstein condensed system with a dark soliton that contains two zero modes corresponding 
to spontaneous breakdown of the U(l) gauge and translational symmetries. Then we evaluate the 
standard deviations of the zero mode operators and see how the mutual interaction between the two 
zero modes affects them. 

PACS numbers: 03.75.Hh, 03.75.Nt, 67.85.-d 


I. INTRODUCTION 

Since the pioneering experiments [1-3], the Bose- 
Einstein condensate (BEC) phenomenon in a trapped ul- 
tracold atomic system has been the central subject of 
many experimental and theoretical studies. We view 
the subject from the standpoint of quantum field the¬ 
ory, which is the most fundamental dynamical law and 
in which the BEC in an atomic system is interpreted as 
a spontaneous breakdown of the U(l) gauge symmetry. 

The concept of spontaneous symmetry breaking (SSB) 
yields many examples of successful descriptions of nature 
using quantum field theory. In an infinite system with a 
spontaneously broken symmetry, the celebrated Nambu- 
Goldstone theorem [4, 5] implies the existence of a zero 
mode, reflecting the original symmetry. The zero mode 
plays a crucial role in creating and retaining the ordered 
state. 

Study of the BEC in a trapped system highlights the 
importance of treating the operators belonging to the 
zero mode sector carefully because the system has a fi¬ 
nite size, and the zero energy state also stands alone as 
a discrete level. The importance is easily overlooked in 
a homogeneous infinite system, for the zero mode sec¬ 
tor is buried in the continuum labeled by the momentum 
index. The practice of formulation, usually called the 
Bogoliubov theory, is to take an unperturbed Hamilto¬ 
nian up to the second power of a field operator and to 
attempt to diagonalize it by expanding the field opera¬ 
tor in the appropriate complete set. However, when the 
complete set includes eigenfunctions belonging to a zero 
eigenvalue, i.e., when the system has zero modes, the un¬ 
perturbed Hamiltonian of the zero mode sector cannot be 
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diagonalized in terms of creation and annihilation oper¬ 
ators but has the form of free particles expressed by the 
quantum mechanical operators P and Q. We refer to this 
as the free zero mode formulation. Although this part of 
the Hamiltonian is simply neglected for a homogeneous 
infinite system, it cannot be neglected for a finite-size 
system. Then, according to Ref. [6], the phase of the or¬ 
der parameter is definite only for a short time because 
Q is interpreted as a phase operator, and its quantum 
fluctuation is given by AQ ~ t for large t. We also point 
out that there is no criterion for specifying a vacuum 
uniquely, as any energy eigenstate of a free particle has 
infinite AQ. Thus, it is concluded that the free zero mode 
formulation for a finite-size system is inconsistent. 

To address this inconsistency, we proposed a new for¬ 
mulation for spontaneous breakdown of the U(l) gauge 
symmetry in Ref. [7]. The key point there is the in¬ 
clusion of higher-than-third powers of P and Q in the 
unperturbed Hamiltonian, which yields their nonlinear 
equations of motion. We therefore call it the interacting 
zero mode formulation. The stationary Schrodinger-like 
equation with the nonlinear unperturbed Hamiltonian of 
P and Q gives bound states rather than the free one, 
and the energy spectrum becomes discrete, so the ground 
state is identified uniquely as the vacuum. Then AQ is 
independent of time, and we have no inconsistency as 
long as the calculated AQ is small. 

The interacting zero mode formulation not only en¬ 
ables us to describe quantum fluctuation of a zero mode 
properly, but also, when two or more symmetries are bro¬ 
ken spontaneously and there are multiple zero modes, in¬ 
troduces interactions among the zero modes naturally. In 
this paper, we focus on interactions among zero modes, 
extending the interacting zero mode formulation for a 
single zero mode [7] to that for multiple ones. After giv¬ 
ing a general formulation, we consider, as an example of 
its application, a system consisting of a dark soliton in a 
homogeneous condensate, where two zero modes coexist 
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corresponding to the spontaneously broken translational 
and U(l) gauge symmetries. This example was studied 
in Ref. [8], in which a nonperturbative treatment of the 
zero modes using an effective Hamiltonian was proposed. 
However, the resultant Hamiltonian represents the free 
motion of each zero mode with no interaction, so our 
approach is essentially different from that in Ref. [8]. 

This paper is organized as follows: In Sect. H, we ex¬ 
tend the interacting zero mode formulation for a single 
zero mode, presented in [7], to a general case of multi¬ 
ple zero modes, comparing it with the corresponding free 
zero mode case. In Sect. HI, the general formulation is 
applied to the originally homogeneous system containing 
a condensate and a dark soliton. Two zero modes ap¬ 
pear and interact with each other. We are particularly 
interested in the quantum fluctuations AQi (where i la¬ 
bels each zero mode), which are affected by interactions 
among the zero modes. Section IV presents a summary 
and conclusion. 


II. FORMULATION OF MULTIPLE ZERO 
MODES IN QUANTUM FIELD THEORY 


In Ref. [7], we proposed the interacting zero mode for¬ 
mulation for a finite-size system with spontaneous break¬ 
down of the U(l) gauge symmetry. The main motivation 
for the formulation is that the quantum fluctuation in 
the phase of the order parameter cannot remain small in 
the conventional free zero mode formulation, whereas its 
smallness is the starting assumption of the formulation. 
In this section, we extend the new formulation to cases 
of multiple zero modes. 

We suppose a Hamiltonian with global U(l) gauge 
symmetry. 


H = 





( 1 ) 


where I4x, w, /i, and g are the external potential, atomic 
mass, chemical potential, and repulsive coupling con¬ 
stant [g > 0), respectively. Throughout this paper, we 
set h to unity. The bosonic field operator obeys the 
canonical commutation relations {x',t)\ = 

6 {x — x'), ['ip{x,t),'ip{x',t)] = 0. When the U(l) gauge 
symmetry is broken spontaneously, ip is divided into co¬ 
herent and incoherent parts as ip = ^+(p, according to the 
criterion (0|(,5|0) = 0. The coherent part, or the order 
parameter is related to the total number of conden¬ 
sates, Nq = J dPx |.Jp, and the vacuum |0) is determined 
self-consistently later. The total Hamiltonian is rewritten 


in terms of ip as H = Hi + H 2 + H 3 + H 4 , where 


Hi= ! (fx -f Hex - M + sICP )? 


H 2 = 


1 


1 




^ Cp + —pAip + -p^A4*p''‘ 


Hs = g 


I A 


2^ ^ 2’ 
+ h.c.. 


H-h.c., (2) 

, ( 3 ) 

( 4 ) 

Hi = ^ jd^x p^’p'^pp , (5) 

where C = —'V'^f2m + I4x — M and M = gp^ . 


A. FREE ZERO MODE FORMULATION 


In the conventional approach, one chooses Hi + H 2 as 
the unperturbed Hamiltonian assuming small p : 

Hii = Hi+H 2 . (6) 

Because the vacuum of this Hamiltonian is time- 
independent, the field division criterion (0|(^(a;,t)|0) = 0 

and the Heisenberg equation yield Hi = 0. This implies 
that ^ should satisfy the Gross-Pitaevskii (GP) equation 
[ 11 ], 

(-|^ + Kx-/i + 5l^p)e = 0. (7) 


To diagonalize H 2 [6, 9], we introduce the Bogoliubov-de 
Gennes (BdG) equation [12, 13] Ty„ = with the 

doublet notations 


T = 





( 8 ) 


Now, we restrict ourselves to cases where all the eigenval¬ 
ues are real, which implies that the system is dynamically 
stable, in order to diagonalize the excited modes. The di- 
agonalization of a system with complex modes (that is, 
a dynamically unstable system) is discussed in Ref. [10]. 
We consider that some symmetries in addition to the 
U(l) gauge symmetry are spontaneously broken and we 
need two or more eigenfunctions belonging to zero eigen¬ 
value, ie., Tyo.i = 0 where yo,i = (/*, -/*)‘, with a 
label i. For the sake of completeness, one has to intro¬ 
duce an adjoint function y-i^i = {hi, h*)* to each ygy, 
which satisfies Tyo.i = and 2 JdPxh* fi = l, with 

normalization constants li. Let us expand p{x,t) by the 
BdG complete set as p{x,t) = pii{x,t) p^y^{x,t), where 


po{x,t) = ^ ^-iQi{t)fi{x) + Pi{t)hi{x) , 

i—z.m. 

pey{x,t) = ^ ^de{t)ui{x) + dl{t)vg{x) , 
£—eyi. 


(9) 

( 10 ) 


where “z.m.” and “ex.” represent summations over the 
zero and excitation modes, respectively. The operators 
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satisfy [Qi{t), Pj{t)] = iSij, [ae{t),al,{t)] = 6u> , and the 
vanishing ones otherwise, where Qi{t) and Pi{t), also 
called the zero mode operators or quantum coordinates, 
are hermitian. Substituting the expansions (9) and (10) 
into Eq. (3), we obtain 


fluctuations cannot be kept small, and the assumption 
of small ifo has to be abandoned. This recognition is the 
starting point of Ref. [7]. Gathering all the terms con¬ 
sisting only of (fo in the total Hamiltonian, we introduce 
the new unperturbed Hamiltonian 


Ho — H 2 — ^ 21 ^ ■ (H) 

i—z.m. ^ £—eK. 

The unperturbed Hamiltonian (11) is represented by the 
sum of the harmonic part of the excitation for each i 
and the free part of the quantum coordinates for each i. 
Because of the latter contribution, we call this conven¬ 
tional formulation the free zero mode one. As a mat¬ 
ter of course, there is no interaction among the zero 
modes in the Hamiltonian (11). The vacuum of the 
total system would be the ground state of Eq. (11): 
|0) = Hi |*^)ex ’ ■''^here |0) . is the ground state in the 
ith zero mode sector satisfying the field division criterion 

,(0|Q,|0)^= ,(0|A|0)^ = 0. 

However, the choice of |0)^ causes problems. First, li 
may be negative, which corresponds to having a system 
with a negative “mass,” and there is no ground state. 
On the other hand, for li > 0, the lowest eigenstate is 
the zero “momentum” state, A|0)^ = 0. As a result of 

the uncertainty relation, the standard deviation of Qi, 

denoted by AQ* = y^j( 0|(5||0). - i(o|(5i|o)^, diverges. 
This divergence immediately conflicts with the starting 
assumption of small ip. It also yields an unphysical sit¬ 
uation in which some physical quantities such as the to¬ 
tal number density also diverge, which will be seen in 
Eq. (29). One way to resolve this contradiction might 
be to choose a wave packet state as the vacuum [6]. Al¬ 
though AQi is finite at t = 0, it is again divergent after 
a long time because of the collapse of the wave packet; 
that is, Qi{t) = (5i(0) + Pi{Q)t and AQi oc t for large t. 
It is argued in Ref. [6] that Q, acting as the phase oper¬ 
ator of the order parameter, yields a divergent AQ{t) as 
t goes to infinity, and that the phase diffuses. All of the 
above pathological properties are rooted in the fact that 
the free Hamiltonian has a continuous spectrum. 

We note that the difficulties concerning quantum fluc¬ 
tuations of zero modes are veiled in the Bogoliubov ap¬ 
proximation, in which the original creation and annihila¬ 
tion operators associated with the eigenfunction belong¬ 
ing to zero eigenvalue are replaced with classical numbers 
in the field ij;, or the zero mode operators in (p are simply 
neglected. 


B. INTERACTING ZERO MODE 
FORMULATION 

Although the choice of the unperturbed Hamilto¬ 
nian (6) is based on the assumption of small tp, or, more 
precisely, small ipo and small (p^^, (0|(5f(t)|0) is diver¬ 
gent, which indicates large (po- Thus, the zero mode 


Hu = H 2 + AH, (12) 


where Hi = 0, for the same reason as in the previous 
subsection. The additional component, AH, is the sum 
of the third and fourth powers of the zero mode operators 
and their counter terms, 


i—z.m.. 


(13) 


The superscript QP indicates that all the terms consist¬ 
ing only of Qi and Pi are picked up. We set up the 
stationary Schrodinger-like equation, 

(i. = 0,l,2,---). (14) 

As was seen in Ref. [7] for the model with the single 
zero mode and will be seen in the example with two zero 
modes in the next section, the eigenequation (14) is a 
type of bound state problem and yields a discrete spec¬ 
trum, in contrast to the free zero mode formulation in 
the previous subsection. It is quite natural to take the 
whole unperturbed vacuum |0) = l^to) ® where 

[tko) is the ground state of the above equation. The un¬ 
known parameters dpi and 5vi involved in H^^ should 
be determined so as to satisfy the field division criterion 

(4'o|ft|4'o) = (4'o|A|4'o) = 0 (15) 


in a manner consistent with Eq. (14). 


Substituting the expansion (9) into Eq. (12), we gather 
it as Hu = where 


= - SviQi, 


hQP-. 

^u ,2 — 


^.^r = 2Re 


-1 ’ 


(16) 

(17) 


'^■^9jkiQjQkQEpP9jkiQj\.Qki 
Hi^jOlPtQkQj P9jjPfiPi 

— ^Cg^i^gPk{Qj,Pi}+DejkePjPkPe , (18) 


tjQP- 

-“«,4 — 


A, 


ijki 


QiQjQkQ^ Im[H2j7i;^] QiQj\^Qk^ Pi} 


— Re[C'ijfcf] QiQjPkPi- 


C'. 


Ijki 


{Q^.Pk}{QJ,Pi} 


-Im[A,fed {4, PpPkPi+^PPjPkPe , 


( 19 ) 
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with 


II 

/ dxjirghfi. 

^ijki — 9 

1 dxf:f;hhg, 

II 

J dxf^*hkhi, 

^ijki = 9 

j dxf*fjhlhi, 

II 

j dxf*h*hkhe, 

^ijki — 9 

J dxh*h*hkhg. 


( 20 ) 


In Eqs. (18) and (19), we define {0i,02} = O 1 O 2 + 
O 2 O 1 , and the dummy indices should be summed over. 
A remarkable consequence of the present formulation is 
that we have cross-terms among the different zero mode 
operators in the unperturbed Hamiltonian, namely, in¬ 
teractions and mixings among them. In other words, 
the Hamiltonian is an effective Hamiltonian in the 
zero mode sector governing the dynamics of the conden¬ 
sate and is uniquely derived from the original Hamilto¬ 
nian (1). 



FIG. 1. (Color online) Plots of the zero modes and their adjoint 
eigenfunctions. Red bold and thin solid lines indicate fg and hg, 
respectively. Green bold and thin broken lines indicate fx and hx, 
respectively. 


III. APPLICATION TO HOMOGENEOUS 
SYSTEM WITH DARK SOLITON 

In this section, we consider a system consisting of a 
dark soliton in a homogeneons system, which is described 
by the Hamiltonian in Eq. (1) with Vex = 0 and has two 
zero modes. The GP solution of a one-dimensional dark 
soliton is 


(22) and (23), respectively, may be interpreted as the 
phase and position operators of the soliton. 


i’ix) = - iQe^{.x) + Qx 

- + Qx)e-^^o. 


di{x) 

dx 


(27) 


^(x) = vfiT^tanh{K(a; - xo)} , = gno , (21) 


where n = y/mgriQ, and ng is the bulk density of the 
condensate. Hereafter, we set xq = 0 and no = 1 for 
the sake of simplicity. The zero modes and their adjoint 
eigenfunctions are 


fe = ^{x), 


fx = , 

ax 

hg = [tanh(«;x)-|-Kx {l —tanh^(nx)}] , 

ZIj 


h 


X 



h = 


5 

L ’ 


Ix 


4k ’ 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


where the subscripts 9 and x denote the 17(1) gange and 
translational modes, respectively. We plot the four func¬ 
tions for later convenience in Fig. 1. Note that although 
in Ref. [8], which considered a nonperturbative treat¬ 
ment of the two zero modes, the system is put in an 
artificial box with a size L and the boundary condition 
fg{—L/2) = fe{L/2) = 0, onr boundary condition is an- 
tiperiodic, fg{—L/2) = —fg{L/2) ^ 0. The operators 
Qe and Qx^ which are associated with the eigenfunctions 


However, this interpretation is true only when Qg and 
Qx are small. Our approach does not rely on it. 


A. FREE ZERO MODE APPROACH FOR DARK 
SOLITON 


In the conventional treatment, the unperturbed Hamil¬ 
tonian (11) for this system is described by 


Hg = 


p2 

xg 


212 


p2 

X 

2/7 


ujialdg , 


(28) 


with positive Ig and negative Ix- As mentioned above, 
the ground state for the U(l) gauge zero mode induces 
phase diffusion, and there is no ground state for the 
translational zero mode because of the negative “mass,” 
Ix < 0. In Ref. [8], Eq. (28) is derived from the classical 
Lagrangian for the collective coordinates as the starting 
point of the nonperturbative treatment, and the U(l) 
gauge zero mode sector is restricted to a subspace with 
a definite total number of atoms, assuming a large phase 
fluctuation. After this procedure, the calculated AQx 
with a Gaussian wave packet state grows with t. As can 
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be seen from the total number density 


(^ 0 \ip^{x, t) |0) 

= [l + AQl{t)] \ax)\^ + AQlit) 


d^jx) 2 _ 

dx 

(29) 


AQx is related to the number of atoms that fill the center 
of the soliton. It has been indicated that AQx{t) grows 
with t, which is quantum depletion. A similar result of 
quantum depletion for a system with an attractive inter¬ 
action having a bright soliton was obtained in Ref. [14], 
although then the “mass” Ix is positive. 

The difficulty of the free zero mode approach lies in the 
the negative “mass” problem, the absence of a ground 
state for the unperturbed Hamiltonian (28), and the fact 
that AQx(t) increases to infinity with time. 


B. INTERACTING ZERO MODE APPROACH 
FOR DARK SOLITON 

In the interacting zero mode approach, we can ob¬ 
tain the natural vacuum, which causes neither phase dis¬ 
sipation nor the negative “mass” problem, by solving 
the Schrodinger-like equation (14) with the unperturbed 
Hamiltonian (12). 

First, one has the nonperturbative Hamiltonian con¬ 
sisting only of Qi and Pi with i equal to either 0 or x, 
denoted by {i = 9, x). In addition, there are in¬ 

teractions between {Qg , Pg} and {Qx , Px}, on which 
we focus our attention. We seek the dominant contri¬ 
bution of the interaction terms in the present model. 
Equations (22)-(25) and Fig. 1 show that fg and hg are 
odd functions with respect to the variable a:, whereas 
fx and hx are even ones. The indices of the 

non-vanishing cross-term coefficients in Eq. (20) must be 
(6,9, X, x) in random order. We consider the limit of large 
L, much larger than the coherent length. Then the mag¬ 
nitude of hx is small, i.e., 1/L, so the D and E terms 
can be neglected. The function fx peaks sharply around 
a; = 0, and the contributions of the A and B terms are 
small. Thus, the dominant contributions come only from 
the C and C terms. One can neglect the third-power 
terms with respect to the zero mode operators because 
their contributions are small compared with those of the 
fourth-power terms owing to the field division criterion. 
Consequently, we have the approximate Hamiltonian 

+ H2^Px+^\Cggxx\QlP^x ■ ( 30 ) 

Here the last term represents the dominant interac¬ 
tion between the U(l) gauge and translational zero 
modes. Solving the Schrodinger-like equation (14) with 
the Hamiltonian (30) numerically, we find the ground 
state or vacuum in the zero mode sector. To illustrate 
the significance of the mutual interaction, namely, the 
last term in Eq. (30), we plot the ground state distribu¬ 
tion |4'o((3e, Qx = 0)P with and without it in Fig. 2; the 


difference between the plots is striking. The sharp distri¬ 
bution in the presence of the mutual interaction can be 
understood from the fact that the last term, behaving as 
Q'g, serves as an additional harmonic potential for Qg. 



- 4 - 3 - 2-101234 



FIG. 2. (Color online) Ground state distribution I’J'olQsiQi = 
0)p with and without the mutual interaction for g = 1 and 
L = 1000. (Red) bold solid and thin broken lines indicate the 
distributions for and + hS,x , respectively. Axes are 

scaled by ig = Igl2AAgegg (see Ref. [7]). 

One can evaluate the standard deviations AQi using 
the ground state distribution obtained above. The re¬ 
sults versus the parameter g are presented in Fig. 3, 
which shows that AQg and AQx decrease according to 
the power law g~°‘% with ag = 0.092-•• and ax = 
— 0.192-It is natural that AQx decreases as g in¬ 
creases, as the coherent length is proportional to g~^^^ 
and the soliton becomes sharper for larger g. We empha¬ 
size that the mutual interaction between the two zero 
modes qualitatively changes the g dependence of AQg: 
There would be no g dependence without the mutual 
interaction, as in a normal condensate [7]. The suppres¬ 
sion of AQg because of the mutual interaction may be 
rephrased as the condensate phase becoming rigid be¬ 
cause of the presence of a soliton with a tt phase kink. 

So far we have considered a homogeneous limit, he., 
the large L limit. We turn our attention to finite-L sys¬ 
tems and estimate the L dependence, for L is a control¬ 
lable parameter in experiments, although in that case 
the translational symmetry is broken not spontaneously 
but explicitly. The calculations of the interacting zero 
mode formulation presented above are straightforwardly 
applied to a cylindrical system with a circumferential 
length 2L and two dark solitons, or to a system con¬ 
sisting of a dark soliton confined in a box of size L. For 
the latter, although our antiperiodic boundary condition 
is not realistic, we expect that the calculations reflect the 
effects of finite L qualitatively. The L dependence of AQi 
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FIG. 3. (Color online) g dependences of standard deviations AQg 
and AQa; with m = 1 and Nq = 10^. (Red) bold and (green) thin 
solid lines indicate AQg and AQj; for the ground state of Hu^ 
in Eq. (30), respectively. Broken and dotted lines indicate the 
same quantities for the ground state of the Hamiltonian + 

hS,x , namely, for the ground state without the mutual interaction 
between the two zero modes. 

is shown in Fig. 4, which shows that AQg decreases as 
with f3g = —0.440 • • •, whereas AQx increases as 
with Px = 0.115---. It is not surprising that AQg de- 



FIG. 4. (Color online) Effect of system size on standard deviations 
AQg and AQx with m = 1 and g = 1. Each line indicates the same 
quantity as in Fig. 3. 

pends on L, as the eigenfunction fe{x) in Eq. (22) extends 
over the entire range from —Lj^ to Lj^. On the other 
hand, the fact that AQx also clearly depends on L is 
puzzling at first glance, because the eigenfunction fxix) 
in Eq. (23), with which the position operator of the soli- 


ton Qx is associated, has a sharp distribution around the 
center of the soliton and is not affected by L. The puzzle 
can be resolved as follows: Whereas Qx describes the po¬ 
sition of the soliton, its conjugate partner Px corresponds 
to its momentum or velocity and is size-dependent, as the 
spread of hx{x) is of order L [see Eqs. (9) and (25)]. It 
is well known that although a standing soliton has a tt 
phase kink, a moving one has a smaller kink [15]. In view 
of this, the global character of hx{x), as seen in Eq. (25) 
to untwist the phase kink, could be understood. Because 
Qx and Px are canonically conjugate to each other, and 
one has the uncertainty relation AQx ■ APx ~ 1/2, the 
standard deviation AQx depends on L as a consequence 
of the L dependence of APx. The size effects may be 
observed in experiments on finite-size systems in which 
L is controlled. 


IV. SUMMARY 

Considering that the zero mode is the essence of SSB 
and that its quantum fluctuation must be treated prop¬ 
erly, we adopted the interacting zero mode formulation 
and extended it from a single zero mode system to a 
general case of multiple zero modes. It yields the effec¬ 
tive Hamiltonian of a pair of canonical operators for each 
zero mode, the spectrum of which is discrete, as in the 
case of a single zero mode system, and introduces inter¬ 
actions among the zero modes naturally and definitely. 
The physical picture of zero modes interacting with each 
other is quite new. 

As an application of the new formulation, a system of 
size L with a dark soliton is considered. In the large L 
or homogeneous limit, there are two zero modes corre¬ 
sponding to spontaneous breakdown of the U(l) gauge 
and translational symmetries. We investigated this sys¬ 
tem by performing calculations. The vacuum is obtained 
uniquely, and the standard deviations for the zero mode 
operators AQi can be evaluated. The mutual interaction 
between the two zero modes influences the ground state 
distribution and therefore AQg, and its effect is seen in 
the way that AQi depends on the coupling constant g. 

As the trapped ultracold atomic system has a finite size 
L in real experimental situations, we also studied the L 
dependence of AQi, keeping L finite in our calculations. 
The results may be checked in experiments on a cylin¬ 
drical system with two dark solitons or a dark soliton 
system confined in a finite region. 

As a characteristic of soliton physics, the behavior of 
the translational zero mode of the dark soliton is ex¬ 
pected to correlate with its velocity. A study of the cor¬ 
relation is a future work. 
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